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Van Hove Singularities in the Quark-Gluon Plasma 

Markus H. Thoma^* 

'^Theory Division, CERN, CH-1211 Geneva 23, Switzerland 

General arguments as well as different approximations for the in-medium quark propagator in a quark-gluon 
plasma lead to quark dispersion relations that exhibit a minimum in one branch (plasmino) . This minimum causes 
Van Hove singularities in the dilepton production rate and mesonic correlators, which might have observable 
consequences. 



1. Introduction 

In 1953 Van Hove discussed singularities in 
the density of states, so-called Van Hove singular- 
ities, in solid state physics. The density of states 
of a system, given by 



6{ljJ - UJn{k)) 



(1) 



with the energy eigenstates a;„(fe), can be ex- 
pressed by the surface integral ||] 



(27r)3 |Vw„(fc)r 



(2) 



Here the quantity in the denominator |Vaj„(fe)| 
can be identified with the group velocity. Due to 
symmetries in a crystal the group velocity van- 
ishes at certain momenta, resulting in a divergent 
integrand in (H). This divergence is integrable in 
3 dimensions, leading to a finite density of states. 
In lower dimensions, however. Van Hove singular- 
ities appear. For example, a 2-dimcnsional elec- 
tron gas shows logarithmic singularities, which 
have been discussed in connection with high- Tc 
superconductors [||. 

Here we want to discuss the role of Van Hove 
singularities in a quark-gluon plasma. We will 
argue that the in-medium quark dispersion rela- 
tion consists of two branches of which one has 
a minimum at some finite value of the momen- 
tum. This leads to a vanishing group velocity for 
the collective quark modes. Interesting quantities 
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such as the production rate of low mass lepton 
pairs and mesonic correlators depend inversely on 
this group velocity. Therefore these quantities, 
which follow from self energy diagrams containing 
a quark loop, arc affected by Van Hove singular- 
ities, which might have observable consequences. 

In the next two sections we will show the ori- 
gin of Van Hove singularities using an effective 
quark propagator for the quark loop in the self 
energies under consideration. First we will use 
the so-called hard thermal loop (HTL) resummed 
quark propagator. After that we discuss the im- 
plications of a quark propagator considering the 
presence of a non-vanishing gluon condensate in 
the deconfined phase as indicated by lattice QCD 
calculations. In section 4 we will argue that the 
minimum in the quark dispersion relation and 
thus the appearance of Van Hove singularities is 
a general property of massless fermions at finite 
temperature. Finally we will discuss the chances 
and problems to observe Van Hove singularities 
and to prove in this way the presence of decon- 
fined collective quarks in relativistic heavy ion 
collisions. 

2. Hard thermal loop approximation 

The HTL resummation technique has been de- 
veloped by Braaten and Pisarski to cure some se- 
rious problems of perturbation theory for gauge 
theories at finite temperature jij . Using only bare 
propagators and vertices naive perturbation the- 
ory lead to gauge dependent and infrared diver- 
gent results for physical quantities. A famous ex- 
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ample is the damping rate of a long-wave, col- 
lective gluon mode in the QGP, which was found 
to depend on the gauge choice and became even 
negative corresponding to a plasma instability in 
covariant gauges jsj. Braaten and Pisarski real- 
ized that this undesirable behaviour is due to the 
fact that higher order loop diagrams contribute 
to lower order in the coupling constant at finite 
temperature. Distinguishing between a soft mo- 
mentum scale gT, where g is the gauge coupling, 
and a hard scale of order T, Braaten and Pis- 
arski isolated the dangerous diagrams. They are 
given by one-loop self energies and vertices, which 
contain only hard loop momenta and are related 
to each other by Ward identities. For example, 
in the case of the polarization tensor in QED or 
QCD a result proportional to g^T^ was obtained^. 
Resumming these self energies and vertices within 
the Dyson-Schwinger equation effective propaga- 
tors and vertices can be constructed. These effec- 
tive Green functions have to be used if all legs of 
the Green function under consideration are soft. 
Otherwise bare propagators and vertices are suffi- 
cient. Physical quantities calculated by applying 
this HTL improved perturbation theory, such as 
the gluon damping rate turn out to be gauge 
invariant. At the same time medium effects such 
as Debye screening, leading to an improved in- 
frared behaviour, and Landau damping are in- 
cluded. For a review of the HTL resummation 
method and its applications to the physics of the 
QGP see Ref.|]. 

Let us now consider the HTL approximation for 
the quark propagator in the QGP. Assuming that 
the temperature is much larger than the mass, 
which holds at least for up and down quarks, we 
can neglect the quark masses. Then the most 
general expression for the fermion self energy in 
the heat bath can be written as [fol 



(3) 



using the notation K = {ko,k), k = \k\. The 
scalar quantities a and b are functions of the en- 
ergy fco and magnitude k of the three momentum. 



^This result has been found before in the high tempera- 
ture approximation jij and much earUer aheady using the 
semiclassical Vlasov equation 



They are given by traces over the self energy 
1 



fcoM7oS)], (4) 
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b{ko, k) ^ ^ [K^ tr(7o S) - ko tr{K^^^ S)] , 
which read in the HTL approximation 



tr(7oS) = 2™: 
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(5) 



with the effective thermal quark mass = 
g^T^ /&. Note that the general ansatz (^) is chi- 
rally invariant in spite of the occurrence of an 
effective quark mass [ p^ . Furthermore the quark 
self energy has a non-vanishing imaginary part 
below the light cone (fcg < k"^) which can be re- 
lated to Landau damping for spacelike quark mo- 
menta. 

For massless quarks it is convenient to decom- 
pose the quark propagator into its hclicity eigen- 
states (fc = fc/fc) 0| 



S{K) 



70 - • 7 , 70 + fc • 7 



2D+{K) 2D_{K) 
The zeros w±(fc) of 
D± = (1 + a) i-ko±k)-b, 



(6) 



(7) 



describe the dispersion relation of the particle ex- 
citation with energy a;+ and of a mode q_, 
called plasmino with energy w_ and nega- 
tive ratio of chirality to helicity. The latter is a 
consequence of the medium, breaking the Lorentz 
invariance of the vacuum. 

The HTL quark dispersion relations are shown 
in Fig.l. We observe that both branches start at 
zero momentum at the same energy ll!± (0) = to, 
and approach the bare dispersion oj = k at large 
momenta. The plasmino branch, which is absent 
in the vacuum, has a minimum at k = 0.408 m^. 
The spectral strength of the plasmino decreases 
exponentially for large momenta, indicating the 
purely collective nature of this mode. 

The emission of thermal lepton pairs from the 
QGP proceeds via the decay of a virtual photon 
produced to lowest order by quark-antiquark an- 
nihilation. For low invariant masses of the or- 
der gT medium effects from the collective quark 
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Figure 1. Quark dispersion relation in HTL ap- 
proximation. Also shown are the free dispersion 
relations of massless (dotted line) and massive 
(dashed line) quarks. 



modes have to be taken into account. The annihi- 
lation amplitude is related to the imaginary part 
of the photon self energy containing a quark loop 
by cutting rules According to the rules of 

the HTL method one has to use effective propaga- 
tors and vertices as shown in Fig. 2, if the photon 
energy and momentum are soft. The imaginary 
part of this diagram comes either from the pole 
of the effective quark propagator, corresponding 
to the dispersion relations of Fig.l, or from the 
the imaginary part of the HTL quark self energy 
contained in the resummed propagator (cut con- 
tribution). Therefore we have pole-pole, pole-cut, 
and cut-cut contributions to the dilepton produc- 
tion rate which is proportional to the imaginary 
part of the polarization tensor according to [Of 



dN 



1 ImH^^(P) 



M2 



(8) 



with the QED fine-structure constant a and the 
invariant photon mass = P^. The pole-cut 
and cut-cut contributions involving external glu- 
ons, as can be seen by cutting the HTL quark 




Figure 2. Photon self energy containing effective 
quark propagators and quark-photon vertices. 



self energy, lead to a smooth contribution to the 
The pole-pole term, on the 



dilepton rate |11|. 
other hand, will cause sharp structures as we will 
discuss below. 

Besides the dilepton production rates also tem- 
poral correlators of mesons follow from the dia- 
gram of Fig. 2 These correlation functions 
of meson currents at the same space coordinate 
but different points in Euclidean time have been 
measured on the lattice in the deconfined phase 
(see e.g. Ref.||l^). A clear deviation from the 
free correlator, which is determined using only 
bare quark propagators and vertices in Fig. 2 [ |l6t , 
has been found in particular in the pion chan- 
nel This indicates the importance of quark 
interactions and medium effects in the QGP at 
temperatures close to the critical. Using HTL 
propagators and vertices as in Fig. 2, medium ef- 
fects, namely effective quark masses and Landau 
damping, are taken into account. It is interesting 
to study to what extend the HTL medium effects 
already explain the lattice results. 

The spectral functions of the temporal corre- 
lators are proportional to the imaginary part of 
the quark loop diagram. The only difference com- 
pared to the photon self energy lies in the vertices 
describing the various mesonic channels. For ex- 
ample, in the pion channel the bare vertex is pro- 
portional to 75. In the vector channel, on the 
other hand, the bare vertex is proportional to 7^ 
as in the dilepton case corresponding to Vector 
Meson Dominance. 

In contrast to photons and vector mesons it is 
sufficient to consider a bare meson-quark vertex 
in Fig. 2 for pseudo-scalar mesons. This is due 
to the fact that HTL vertices in the case of the 
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Yukawa theory lead only to higher order contri- 
butions and can be neglected therefore [0 . Since 
HTL vertices are complicated functions of the en- 
ergy and momentum, we will discuss here only 
the spectral function of the pseudo-scalar tempo- 
ral correlator. Actually the consideration of the 
effective vertex does not change the position of 
the Van Hove singularities, which are determined 
only by the minimum in the plasmino dispersion. 
Furthermore in the dilepton rate from the QGP 
the inclusion or neglect of the effective vertex does 
also not alter the magnitude of the rate signifi- 
cantly. This is in contrast to the dilepton rate 
from pion annihilation, where the consideration 
of the Ward identity reduced similar singularities 
in the rate strongly ||l8| . 

The spectral function of the pseudo-scalar cor- 
relator is given by 



(9) 



aps{uj,p) = - Im Xps{i^,p) 

TT 

where the correlation function in momentum 
space contains the quark loop according to 



(27r)3 



Tr 



755'HTL(fco,fe) 



X 75 5'htl(w - fco,p- fe) 



(10) 



Here TVc = 3 is the number of colors, and the 
sum over n denotes the sum over the fermionic 
Matsubara frequencies ko = (2n + l)iTrT. 

The HTL resummed quark propagator is con- 
veniently expressed by its spectral function phtl 

SFiko,k)^-i^oko~7-k) [' dre'=«^ 

Jo 

/oo 
duj pnTLiuJ, k) [1 - npiuj)] e-^\ (11) 
'OO 

where ni?(cj) = 1/[1 -I- exp(w/r)] and 
PHTL(fco-fc) ^ p+(fco,A;)(7o - « fc • 7) 



+ -p_(fco,fc)(7o + «fe-7) (12) 



with 11 



fc^ - k^ 

P±{ko,k) = -^—-^[6{ko-LU±) + 6{ko + LU^)] 



2m^ 



+ (3±iko,k)e{k^-kl). 



(13) 



Here the first part of (|13|) corresponds to the pole 
contribution of the HTL propagator. The second 
part, corresponding to the cut contribution from 
the imaginary part of the HTL quark self energy, 
is given by 



(3±{ko, k) 



(±fco - k) 



M-fco±fc) + mM±l_±^ln^ + ^° 



2k 



k — ko 



+ 



(14) 



Combining to (14) the spectral function of 
the pseudo-scalar correlator can be written as 

ap,(.;,p)=2iV,(e-/^-l) J 

/OO 
dxdx'nF{x)np{x')6{u) — x — x') (15) 
-oo 

x|(l - q • k)[p+{x,k)p+{x' ,q) + p^{x,k)p^{x,q)\ 

+ {l + q-k) [p+{x, k)p_ (x', q) + p- (x, k)p+{x, q)] | , 

where q = p — k. Considering only the pole-pole 
contribution originating from the first equation of 
( [l^ ) all the integrals in (|5|) can be done exactly. 
One finds 



aPP(c.,p = 0) = 



iVc (e"/^ - 1) 



27r2 



m„ 



nj,{uj+(ki)) 



{Lul{k,)-klfkl 



2kV(fci)l 

2 

2Y,nF{u^+{kl))[l-nF{uJ-{kl))\ 

4=1 

{iol{k^)-{k\f) H{kl)-{k\f){kif 

W+{kl) - L0'4k\)\ 



(16) 



-^n|,(^_(fc^)) 



i=l 



{u:l{kl)-(klYf(klf 
2\Lo'_[kl)\ 



Here u>±{k) denote, as before, the quark disper- 
sion relations for the ordinary quark (-I-) and 
the plasmino (-) branch, fci is the solution of 
u) — 2Lo+{ki) = 0, k2 and k^ are the solutions 
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of LO — LO+m) + uj-m) — and uj — 2L0-{kl) — 0, 
respectively. Note that for small momenta the 
last two equations can each have two solutions. 
Furthermore, (^'±{k) = {duj±{x)/dx)\x=k- 

Eq. (16) is shown in Fig. 3 for niq/T = 1, 



I.e., g = v6. The first part of (|16|) corresponds 
to the annihilation of collective quark-antiquarks 
and sets in at the threshold w > 2mg. For large 
energies this is the dominating contribution which 
approaches the result obtained from a bare quark 
propagator (crosses) for uj » m^. The second 
part corresponds to a transition from the upper 
quark branch to the lower plasmino branch. It 
starts at Lu ^ and terminates with an Van 
Hove singularity at w = 0.47 , where the dif- 
ference u!+{k2) — ^-(fcj) has a maximum due to 
the minimum in the plasmino dispersion. The 
third part describes the annihilation of plasminos 
and antiplasminos which starts at the threshold 
u! = 1.86 mg, where another Van Hove singularity 
occurs due the minimum of the plasmino branch 
(^-(fcg). For larger momenta this contribution 
vanishes quickly due to the exponentially sup- 
pressed spectral strength of the plasmino branch. 
The singularities in (^6|) can be integrated lead- 
ing to finite results for the temporal correlator 
14 1 . Also in Fig. 3 the smooth pole-cut and cut- 
cut contributions are shown. In particular the 
latter is suppressed by about an order of magni- 
tude. 

In the case of the spectral function of the vector 
correlator Q and the dilepton production rate at 
zero momentum ]ll[ , for which an effective HTL 
vertex has to be used, the pole-pole and pole-cut 
contributions are similar (apart from some pre- 
factors). In particular the Van Hove singularities 
are at the same positions. The cut-cut contribu- 
tion, however, diverges for small uj covering up 
the low-energy Van Hove peak and leading to an 
infrared divergent vector correlator. 

The temporal correlator follows from its spec- 
tral function by an energy integration (l4| 



Gpsir) 



duj (7pg{uj, p = 0) 



cosh(w(T-/3/2)) 



sinh(c<j/3/2) ' ^^^^ 
where (3 = 1/T and r is restricted to the 
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Figure 3. Pole-pole, pole-cut, and cut-cut contri- 
butions (solid lines) to the spectral function of the 
pseudo-scalar temporal correlator at rUq/T = 1. 
The crosses indicate the free spectral function. 



Euclidean time interval [0,/?]. Using (0) to- 
gether with ( [l6| ) and the corresponding expres- 
sions for the pole-cut and the cut-cut contribu- 
tions the temporal pseudo-scalar and vector cor- 
relators turn out be similar to the free correla- 
tion functions containing only bare quarks. This 
is caused by a competing effect of the pole-pole 
contribution and the pole-cut and cut-cut con- 
tributions. Whereas the pole-pole contribution 
reduces the spectral function due to the pres- 
ence of the effective quark mass compared to the 
free spectral function, there is an enhancement 
due to the pole-cut and cut-cut contributions de- 
scribing higher order diagrams involving external 
gluons. Surprisingly the two effects compensate 
each other in the temporal correlator almost com- 
pletely. Hence the deviation from the free correla- 
tor, observed on the lattice, cannot be attributed 
to HTL medium effects 

3. Gluon condensate quark propagator 

The problem connected with the above calcu- 
lations relies in its perturbative nature. Since 
the coupling constant g is not small in realistic 
situations, the extrapolation to relativistic heavy 
ion collisions is questionable. Furthermore due 
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to collinear and magnetic divergences the conver- 
gence even of the HTL improved perturbative se- 
ries at finite temperature is upset at least for low 
mass M ~ g^T dileptons 

In Ref. another effective quark propagator 
has been constructed by taking into account the 
gluon condensate measured in lattice QCD cal- 
culations above the phase transition. In this way 
non-perturbative effects are included, which allow 
to study realistic temperatures addressed on the 
lattice. For this purpose we calculate the quark 
self energy from a one-loop diagram using a non- 
perturbative gluon propagator that contains the 
gluon condensate analogously to the zero tem- 
perature case |21 . Of course, the approach is 
a purely phenomenological combination of Green 
functions with lattice results. The effective quark 
propagator can then be written as in (||) and (^) 
with 
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a = — 



ill. 
6 

1 



ok' - -ki {£ 



— k k, 
5 



/ {B')t 



, 4 2 fco 



(18) 



where the in-medium chromoelectric, ^)t, and 
chromomagnetic condensates, {B'^)t, are taken 
from lattice calculations p^ ]. 

In Fig. 4 the dispersion relations following from 
the poles of this propagators are shown at T = 
2Tc |2^]. The important point here is that, al- 
though the HTL resummed and the gluon con- 
densate quark propagators are completely differ- 
ent, the dispersion relations show the very same 
behaviour. In particular we observe again the 
minimum in the plasmino branch. Therefore we 
expect again Van Hove singularities in the dilep- 
ton rate following from this quark propagator. 
In a first step we neglected an effective quark- 
photon vertex, related to the gluon condensate 
quark propagator by Ward identities How- 
ever, as we discussed above, this will not change 
the existence and positions of Van Hove singulari- 
ties following solely from the quark dispersion re- 
lations. In Fig. 5 the dilepton production rate at 
zero photon momentum versus the invariant pho- 
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Figure 4. Quark dispersion relations according to 
the gluon condensate quark propagator at T = 

2Tr. 



ton mass M from this investigation is shown 
for different temperatures. The interpretation of 
the various curves is identical with the one in the 
last section. However, in this case no pole-cut 
and cut-cut contribution are present as the quark 
self energy containing the gluon condensate has 
no imaginary part [20|] . 

4. General quark dispersion relations 

The two completely different approximations 
for the quark propagator, discussed in the last two 
sections, lead to the same qualitative behaviour of 
the dispersion relations. Hence one may speculate 
that the dispersion relations shown in Fig.l and 
Fig. 4 are a consequence of the most general full 
propagator for massless fermions at finite tem- 
perature. Analyzing the most general expression 
for the in-medium fermion propagator (^) and (0) 
one can show the following general properties of 
the dispersion relations (1) there are always 
two branches; (2) both branches start at the same 
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Figure 5. Dilepton production rate at zero photon 
momentum versus the invariant photon mass M 
foUowing from the gluon condensate quark prop- 
agator. 

energy, i.e. effective fermion mass, at zero mo- 
mentum; (3) the both branches have an oppo- 
site slope at zero momentum; (4) the branch cor- 
responding to the usual fermion approaches the 
free dispersion relation, oj = k, for large momenta 
k. Assuming that the plasmino branch also ap- 
proaches the free dispersion at large momenta, 
as it does in the two independent examples dis- 
cussed above, the plasmino branch must always 
possess a minimum. So far this assumption has 
not been proven. Maybe future lattice calcula- 
tions might be capable of investigating the quark 
dispersion relations p6| . However, it appears to 
be reasonable that the plasmino minimum is a 
general feature leading to Van Hove singularities 
in the dilepton production rate and mesonic tem- 
poral correlators. 

5. Consequences for relativistic heavy ion 
collisions 

Now the important question is whether the Van 
Hove singularities in the QGP can be observed 



in relativistic heavy ion collisions. One possibil- 
ity would be to look at the low-mass dilepton 
spectrum, which has been investigated already 
at SPS energies |2^. Here no sign for interest- 
ing structures due to Van Hove singularities has 
been found. However, hydrodynamic simulations 
IH] show that the QGP contribution, if there is 
a QGP phase at all as indicated at SPS are 
suppressed at least by an order of magnitude com- 
pared to the hadronic contributions. This situa- 
tion will change at RHIC and LHC, where the 
QGP contribution should dominate the thermal 
emission from the fireball. If low-mass dileptons 
should be investigated at RHIC or LHC, interest- 
ing structures coming from the Van Hove peaks 
and the gap in the dilepton rate might appear. 
From a comparison of the equation of state found 
in lattice calculations and quasiparticle models 
|3C| ] effective quark masses of the order 500 MeV 
are expected. Therefore the structures may show 
up at invariant masses below about 1 GeV. 

To what extent these distinct structures will 
survive in heavy ion collisions depends on how 
much of it will be covered and smeared out by 
higher order effects, e.g. bremsstrahlung and 
damping. Also finite momenta ||2j,^ and the 
space-time evolution of the fireball will wash out 
the structures to some degree. These effects are 
an interesting subject for future investigations. 
Anyway, if a non-trivial structure should be ob- 
served in the low-mass dilepton spectrum it would 
provide an unique signal for the QGP formation, 
since one does not expect a similar signal from 
the hadronic phase. Due to in-medium collisions 
between hadrons the structures from the hadronic 
phase such as the p-peak are washed out [Q. 
Hence structures in the low mass dilepton spec- 
trum due to Van Hove singularities would not 
only indicate a deconfined phase, but also the ex- 
istence of collective quark modes in the QGP. 

6. Conclusions 

The quark dispersion relations in the QGP have 
been studied using two different approximations 
for the quark propagator. First the HTL method, 
which is based on the resummation of the pertur- 
bative quark self energy obtained in the high tem- 
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perature limit has been adopted. Secondly the 
gluon condensate measured on the lattice above 
the critical temperature has been included in the 
effective quark propagator. The poles of the effec- 
tive quark propagators, following from these com- 
pletely different approaches, determine quark dis- 
persion relations which show the same behaviour. 
In both approaches two branches, of which one 
shows a minimum, namely the plasmino branch, 
have been found. There are strong indications 
that this is a general feature of massless fermions 
at finite temperature. 

The minimum in the plasmino branch leads to 
a vanishing group velocity of the plasma modes 
and therefore to a diverging density of states in 
the low-mass dilepton production rate and in the 
spectral function of mesonic correlators. There- 
fore Van Hove singularities appear in these quan- 
tities. The dilepton production might serve as a 
promising signature for the QGP formation in rcl- 
ativistic heavy ion collisions. Whether these Van 
Hove singularities can be observed in the dilepton 
spectrum from a dynamical QGP, produced pos- 
sibly in the fireball of a nucleus-nucleus collision, 
is an open question. Anyway it will be worthwhile 
to look for non-trivial structures, indicating the 
existence of deconfined, collective quark modes, 
in the low-mass dilepton spectrum at RHC and 
LHC. 

The meson correlation functions, on the other 
hand, can be compared to lattice calculations, 
which allows to extract informations on the non- 
perturbative nature of the QGP. Future QCD 
lattice calculations might be able to investigate 
directly the quark dispersion relations and the 
spectral functions of the mesonic correlators |33| , 
which contain much more informations than the 
correlators themselves. Summarizing, the verifi- 
cation of Van Hove singularities in the QGP ex- 
perimentally as well as on the lattice will be an 
interesting challenge. 
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